Testing Scalar- Tensor Gravity Using Space Gravitational- Wave 

Interferometers 



Paul D. Scharre and Clifford M. WillB 

McDonnell Center for the Space Sciences, Department of Physics, 
Washington University, St. Louis, Missouri 63130 



Abstract 

We calculate the bounds which could be placed on scalar-tensor theories 
of gravity of the Jordan, Fierz, Brans and Dicke type by measurements of 
gravitational waveforms from neutron stars (NS) spiralling into massive black 
holes (MBH) using LISA, the proposed space laser interferometric observa- 
tory. Such observations may yield significantly more stringent bounds on 
the Brans-Dicke coupling parameter uj than are achievable from solar sys- 
tem or binary pulsar measurements. For NS-MBH inspirals, dipole gravita- 
tional radiation modifies the inspiral and generates an additional contribu- 
tion to the phase evolution of the emitted gravitational waveform. Bounds 
on u) can therefore be found by using the technique of matched filtering. We 
compute the Fisher information matrix for a waveform accurate to second 
post-Newtonian order, including the effect of dipole radiation, filtered using 
a currently modeled noise curve for LISA, and determine the bounds on u) 
for several different NS-MBH canonical systems. For example, observations 
of a IAMq NS inspiralling to a 10^ Mq MBH with a signal-to-noise ratio of 
10 could yield a bound of a; > 240, 000, substantially greater than the current 
experimental bound of w > 3000. 

I. INTRODUCTION AND SUMMARY 



The observation of gravitational waves from inspiralling compact binaries is expected to 
provide an excellent testbed for theories of gravity alternative to general relativity (GR) (see 

for a review). In earlier work [H, we have shown that the measurement of gravitational 
waves from neutron star (NS) - black hole (BH) mergers by ground-based detectors such as 
the Laser Interferometric Gravitational- Wave Observatory (LIGO) could bound the strength 
of a scalar field in scalar-tensor theories of gravity. The bounds estimated in ||^ were 
potentially better than the solar system bounds that were in force at the time (1994), but 
are no longer competitive with bounds that have improved through the use of Very Long 
Baseline Interferometry (VLBI) to measure the deflection of light by the sun. 

In this paper, we analyze potential bounds on scalar-tensor theory made possible by 
plans for space-based gravitational- wave observatories such as the Laser Interferometer Space 
Antenna (LISA), a proposed project of the National Aeronautics and Space Administration 
(NASA) and the European Space Agency (ESA) ||^. We show that the detection of NS 
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inspirals into massive black holes (MBH) by LISA may provide an opportunity to place 
significantly more stringent bounds on scalar-tensor theory than current solar-system or 
binary-pulsar bounds or those achievable by an Earth-based detector. 

Space-based interferometers are being designed to detect gravitational-wave signals in 
a lower frequency band than ground-based detectors such as LIGO, GEO, VIRGO, and 
TAMA, which measure waves from approximately 10 Hz to 10^ Hz. The latter interfer- 
ometers are expected to detect waves from low-mass inspiralling compact binaries (masses 
1 — IOMq). LISA is being designed to detect inspirals and mergers of massive black holes 
(MBH) with masses of 10^ — IO^Mq , among other potential sources Mergers of such 
high-mass objects will produce waves in LISA's proposed sensitive band of frequencies be- 
tween 10"^ Hz and IQ-^ Hz. 

The inspiral of such binary systems is driven by gravitational radiation reaction, which 
imprints onto the phase and amplitude of the emitted gravitational waves information about 
the masses and spins of the orbiting bodies, about the orbital elements and their evolution, 
and about the theory of gravity that is in force. Because broad-band interferometric detec- 
tors are particularly sensitive to the phasing of the waves, accurate determinations of system 
and theoretical parameters are possible 0. 

The simplest version of scalar-tensor gravity is that studied by Jordan, Fierz, Brans and 
Dicke, commonly known as Brans-Dicke (BD) theory |^. In BD, a scalar gravitational field 
(f) is postulated in addition to the spacetime metric g^j,^, with an effective strength inversely 
proportional to a coupling parameter uo. In BD, uo is assumed to be a fixed constant. 
(Generalized scalar-tensor theories that are currently popular byproducts of string theory 
and other unification schemes treat t^; as a function of 0.) As — oo, the effect of the scalar 
field tends to zero and BD — > GR. Current experimental bounds on uj from measurements 
of light bending place u > 3000 (for a review see [^]; for recent light bending measurements, 
see i). 

The details of the behavior of compact binaries in scalar-tensor gravity have been worked 
out by Eardley 0, Will 0, Zaglauer [|1^ and Damour and Esposito-Farese [|lT|. The chief 



difference between general relativity and scalar-tensor gravity is the existence, in the latter, 
of dipole gravitational radiation. In general relativity, for slowly moving systems, the leading 
multipolar contribution to gravitational radiation is quadrupolar, with the result that the 
dominant radiation-reaction effects are at order (v/c)^ relative to Newtonian gravity, where 
V is the orbital velocity. The rate, due to quadrupole radiation in GR, at which a binary 
system loses energy is given by 
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where r, v, and r represent the orbital separation, relative orbital velocity, and radial velocity, 
respectively. We use units in which G = c = 1. The quantities rj and m are the reduced mass 
parameter and total mass, respectively, given hj t] = 17117712/ (mi + ^2)^, and 771 = mi + 7712- 
In BD, this formula is modified by corrections to the coefficients of 0{l/u) (BD also predicts 
monopole radiation, but in binary systems it contributes only to these 0(l/a;) corrections). 

The important modification in BD is the additional energy loss caused by dipole radia- 
tion. By analogy with electrodynamics, dipole radiation is a (v/c)^ effect, potentially much 
stronger than quadrupole radiation. However, in scalar-tensor theories, the gravitational 
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"dipole moment" is governed by the difference si — S2 between the bodies, where Si is a mea- 
sure of the self-gravitational binding energy per unit rest mass of each body. Technically, Sj 
is the "sensitivity" of the total mass of the body to variations in the background value of 
Newton's constant, which, in these theories, is a function of the scalar field: 



'(9(ln mj) 
9(ln G) 



I , (1-2) 
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where G is the effective Newtonian constant at the star and the subscript denotes holding 
baryon number fixed. The energy loss caused by dipole radiation is given by 

(1.3) 

' " \ ^ / 

valid to first order in l/u, where S = si — S2- 

In BD, the sensitivity of a black hole is always sbh = 0.5, while the sensitivity of 
a neutron star varies with the equation of state and mass. For example, s^s ~ 0.12 for a 
IAMq NS with a stiff equation of state, and s^vs ~ 0.19 for a I.IMq NS with a soft equation 
of state. A more detailed discussion of neutron star sensitivities can be found in [lyjl^ in 



the context of Brans-Dicke theory, and in [|TT| in the context of general scalar-tensor gravity. 
For white dwarfs, s\yd < 10~^. 

Ironically, binary black-hole systems are not at all promising for studying dipole gravi- 
tational radiation because sbhi — sbh2 = 0, a consequence of the no-hair theorems for black 
holes. Essentially they radiate away any scalar field, so that a binary black hole system 
in scalar-tensor theory behaves as if GR were valid (see for further discussion). Sim- 
ilarly, binary neutron star systems, such as the Hulse- Taylor binary pulsar PSR 1913-1-16 



and similar systems, are also not effective testing grounds for dipole radiation ||10[. This 
is because neutron star masses tend to cluster around the Chandrasekhar mass of IAMq, 
and the sensitivity of neutron stars is not a strong function of mass for a given equation of 
state. Hence in systems like the binary pulsar, dipole radiation is naturally suppressed by 
symmetry, and the bound achievable cannot compete with those from the solar system [jl3| . 
Thus the most promising systems are mixed: BH-NS, BH-WD, or NS-WD. 

While a leading scientific motivation for LISA is detection of waves from the merger of 
binary massive black holes in the centers of galaxies at cosmological distances, LISA will 
also be able to study waves from the inspiral of low-mass compact objects (black holes, 
neutron stars, white dwarfs) into massive black holes in galaxies and globular clusters. The 
orbits of these inspiralling bodies are expected to be complex: possibly highly eccentric 
because the bodies will have been sent toward the hole by gravitational scattering from 
stars in the surrounding cluster; possibly strongly perturbed by the surrounding stars each 
time the body passes through its apholion; possibly dragged in a precessing orbit by the 
Lens-Thirring effect of a rapidly rotating hole. Analysing the gravitational waves from such 
orbits will be a challenging task. 

However, to study the potential for testing scalar-tensor gravity, we will make a drastic 
simplification of the problem: we will assume that the compact body spirals toward a non- 
rotating hole on an adiabatically shrinking circular orbit, up to the innermost stable circular 
orbit (ISCO) of the black hole, whereupon the body plunges across the horizon on a short 
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timescale. We will assume that the waves from such a quasi-circular orbit can be studied 
for up to a year using LISA. 

With these assumptions, one can make a crude estimate of the bounds on scalar-tensor 
theory that could be possible. The idea is to use matched filtering of a theoretical template 
waveform against the output of the detector. A "match" is a signal whose total accumulated 
phase over the integration time of the experiment matches that of the template within a 
number of radians set by the time resolution of the instrument and the relevant frequency. 
The accumulated phase is given by 



$GW = 2nfdt = f^' 2n{f/f)df , (1.4) 



where / is the gravitational- wave frequency (twice the orbital frequency), and the subscripts 
"i" and "f" denote the values at the beginning and end of the integration. The evolution of 
the frequency is given by 

u = M-\96/5y^/\l + hrf"'u-^/^) , (1.5) 

where Ai = rf/^m is the "chirp" mass, and u = iiAif . The leading term in Eq. ( [1.5| ) is 
the contribution of quadrupole radiation damping from Eq. ( |1 . 1|) . The second term is the 
contribution of dipole radiation damping [Eq. (|1.3|)], where b = (5/48)5^/0;. Integrating 
Eq. ( |1.5|) using only the leading term, we obtain approximate initial and final frequencies 
corresponding to an observation period T leading up to the innermost stable circular orbit 
whose gravitational- wave frequency is /isco = (6^'^^7rmsH)~^; Eq. ( |1.5| ) then yields the 
dipole contribution to the phase given by $D ~ 2.46(5Vtur7)[(l + 16r/T/405m)7/s - 1]. For 
a 1.4 Mq neutron star and a 10^ Mq black hole, the frequency at the ISCO is less than 0.5 
Hz, while, with five million kilometer arm lengths, the minimum time resolution for LISA is 
of order 30 seconds, hence the minimum phase resolution for waves from such a black hole is 
of order 100 radians. Demanding that the accumulated phase offset due to dipole radiation 
reaction be less than this resolution leads to the crude bound on u given by 

4 / 5 \ V 100 \ / T [10^ Mq^^^^ 




->''<'0'[r3) lAsrllTTTl 1—^1 . (1.6) 

We confirm this crude estimate with detailed calculations using matched filtering to 
estimate the parameters of such NS-MBH inspiralling systems, and employing current pro- 
posals for the instrumental noise curve for LISA together with an estimate for confusion 
noise caused by a background of galactic white dwarf binaries. The results are displayed 
in Fig. 1. Shown are the bounds achievable for 1.4M0 neutron stars inspiralling into black 
holes with masses between 1000 and 200,000 Mq, and for integration times of one year, 1/2 
year and 1/4 year leading up to the innermost stable orbit. Notice that the bound decreases 
in almost inverse proportion to black hole mass in agreement with our crude estimate; this is 
a result of the lower intrinsic frequencies of high mass systems, and the consequent accumu- 
lation of fewer cycles in the relevant integration time, leading to less ability to separate out 
the dipole effects from quadrupole phasing. Notice also that, for a given mass, the bound 
scales almost linearly with integration time, again in agreement with our crude estimate. 
Throughout, we assume for concreteness a maximum (and thus) pessimistic value sns = 0.2 
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for the sensitivity of the neutron star, so that S = 0.3. If sns were as low as 0.1, the bounds 
would increase by (0.4/0.3)^ ~ 1.7. For IAMq white dwarf inspirals {swd ^ sbh), the 
bounds would be larger by (0.5/0.3)^ ^ 2.8, provided that Mbh > 2 x 1O*^M0, so that the 
white dwarf reaches the ISCO without tidal disruption. For Mbh < 2 x IO^Mq, integration 
must be cut off as many as a dozen days before the ISCO, with a consequent loss of sensi- 
tivity, so that for a IOOOMq black hole, white dwarf inspiral leads to almost the same bound 
as does NS inspiral. 

The remainder of the paper is devoted to details of the method and calculations. In Sec. 
^ we summarize the basic method of parameter estimation using matched filtering; this 
method is based on foundations laid for gravitational-wave detection by Finn and Chernoff 
and Cutler and Flanagan [^, and applied to specific parameter estimation problems 



for inspiralling systems by Will and Poisson P JT^JTTH . In Sec. |T| we apply the method 
specifically to observations by space interferometers. First we consider LISA. Then we we 
apply the method to a hypothetical follow-on to the LISA mission (which we call SuperLISA), 
whose sensitive band would lie between that of LISA and the ground based interferometers. 
The results are discussed in Sec. |V| 



II. ESTIMATING BINARY SYSTEM PARAMETERS USING MATCHED 

FILTERING 

It is expected that laser interferometric observatories, whether ground-based or space- 
based, will endeavor to detect gravitational waves from binary inspirals using a technique 
known as matched filtering wherein a template consisting of a theoretical gravitational 
waveform is compared to the detector output. The template that is a match with an actual 
signal present with the noise will show a strong correlation with the detector output and 
thereby will "filter" a signal out of background noise. The template is represented by h{t), 
related to the spatial components of the radiative metric perturbation far from the source. In 
the so-called "restricted post-Newtonian approximation" for an inspiral orbit that is quasi- 
circular (that is, circular apart from the adiabatic decrease in separation), we approximate 
h{t) ~ 3?{/io(t)e***^*^}, where /io(^) is a slowly varying wave amplitude; it depends on the 
wave polarization, the location of the source on the sky relative to the detector, the detector 
orientation, and the distance to the source; $(t) is the phasing of the wave, which is a 
function of the evolving orbital frequency; and 3? denotes the real part. The wave amplitude 
/lo(^) is considered only to Newtonian order, i.e. given by the lowest-order, quadrupole 
approximation, because the process of matched-filtering is rather insensitive to changes in 
the amplitude of the wave. Matched filtering is very sensitive to the phasing of the wave, 
however, and thus $(t) is taken to the highest post-Newtonian order possible. The Fourier 
transform of h{t) using the stationary phase approximation is given by 

M/)= n 'U</</max ^2.1) 

[ U , J ^ J max 

where /max is the largest frequency for which the wave can be described by the restricted post- 
Newtonian approximation. For inspiral into black holes, this is often taken to correspond to 
waves emitted at the innermost stable orbit (ISCO) before the bodies plunge toward each 
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other and merge. In the hmit where one mass is much smaller than the other, this frequency 
is given by 

/max = (6=^/Vm)-i. (2.2) 

All the relevant factors related to the polarization of the wave, distance to the source, and 
its position on the sky are included in the amplitude A. After averaging over all angles, we 
obtain 

1 A^5/6 



where is the "luminosity" distance of the source. 

Whereas the amplitude is calculated to Newtonian order, the phase \I/(/) includes higher 
post-Newtonian corrections. To Newtonian order in general relativity, '^{f) is given by 

^U) = ^^ftc -ct^c- vr/4 + ^u-'/' , (2.4) 

where u = irAif and 0c is formally defined as the phase of the wave at the time of coales- 
cence, tc- 

In BD, for quasi-circular orbits, this formula is modified in two ways. (1) To leading 
order, the modifications of quadrupole and monopole radiation can be subsumed in the 
simple replacement of the GR chirp mass A4 in Eqs. ( |2.^ ) and ( |2.4| ) by an effective chirp 
mass given by 

^3/5 

M^-^v'^'m, (2.5) 

where 

Q _ ^ _ ~^ ^2 ~ 2SiS2 

2 + iu ' 

" = ^[^^2(2 + ^) 12(2 + ^) j ' ^^-^^ 

and where Si and S2 are the sensitivities of the two bodies. (2) Via the energy loss in Eq. 
( |1.3| ), dipole radiation introduces another term in the phase "^{f), so that Eq. ( p.4| ) now 
reads 

^{f)BD = 2n ft, - 0, - n/4 + - ^fer^^/s^-^/^) , (2.7) 

where the parameter b is given by 

48 2 + cj ^ ^ 

As uj ^ 00, 6^0. For more details on the gravitational waveform h{t) in BD, see [Q. 
Because we assume a priori that uj > 3000, we shall approximate b ^ (5/48)5^/^7. We 
further add to the phasing formula the PN term and the 1.5 PN "tail" term of general 
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relativity; they serve the useful purpose in the matched filter of breaking the degeneracies 
among the various parameters through their different frequency dependences. While not 
strictly correct in BD, they are individually valid up to corrections of order l/cu. Thus we 
adopt the phasing 
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(2.9) 



where the third and fourth terms inside the square brackets are the PN and 1.5PN terms, 
respectively. Notice that, because m ~ v^, the dipole term is O^l/v"^) compared to the 
quadrupole term, as expected. 

By maximizing the correlation between a template waveform that depends on a set of 
parameters 6"' (for example, the chirp mass Ai) and a measured signal, matched filtering 
provides a natural way to estimate the parameters of the signal and their errors (for discus- 
sion, see 0,0). With a given detector noise spectrum Sn{f) one defines the inner product 



between two signals hi{t) and h2{t) by 

{hi\h2) = 2 — df, (2.10) 

-'o Sn{f) 

where hi{f) and h2{f) are the Fourier transforms of the gravitational waveforms h{t). The 
signal-to-noise ratio (SNR) for a given h is given by 

p[h]^ih\hy/\ (2.11) 

One then defines the "Fisher information matrix" Tab with components given by 

( dh dh\ 



ah 



89" ' de^ 



{2.12) 



An estimate of the rms error, A^°, in measuring the parameter 9°" can then be calculated, 
in the limit of large SNR, by taking the square root of the diagonal elements of the inverse 
of the Fisher matrix, 



Ar = VS«^, S = r-^ (2.13) 
The correlation coefficients between two parameters 9°" and 9^ are given by 



c° = E'^VVS^'^S^^ (2.14) 

In our calculations, the parameters to be estimated will be 0^, fo^o InAI, In 77, and 6, 
where h = hrf/^., and /o is a fiducial frequency characteristic of the detector noise spectrum. 
We follow the method used in 0: combining Eqs. (|2.1| ) and ( ^^ and calculating the partial 
derivatives dh / 89"- for the five listed parameters, we construct the Fisher information matrix 
using Eqs. (|2.1CI| ) and ( p.l2|) . For simplicity, we set h equal to its nominal GR value of zero 



in all information matrix expressions. We then invert the information matrix and evaluate 
the errors in the five parameters, along with the correlation coefficients between A^, rj and 
h. Since the nominal value of h is zero {uj = 00), the error on it translates into a lower bound 
on UJ. Results are calculated for the canonical neutron star mass of I.4M0, and for various 
values of the black hole mass. 
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III. BOUNDS ON SCALAR-TENSOR GRAVITY FROM SPACE-BASED 

INTERFEROMETERS 



A. LISA-type interferometer 

We consider space-based interferometers of the proposed LISA type, with a sensitive 
bandwidth between 10~^ and 10~^ Hz, typical integration times up to one year, and an 
expected noise curve which can be expressed in terms of an overall amplitude 5*0, and a 
function of the ratio x = f / f^: 

Snif) = Sogix) . (3.1) 
Including the LISA instrumental noise curve and an estimate of "confusion noise" from a 



population of galactic white-dwarf binaries ||18|JT9|1 , we adopt a noise curve given by 

So = 4.2 X 10"^^ Hz"^ , 
/o = 10-3Hz, 

g{x) = VlOx-'^/' + 1 + xVlOOO + 3i3.5a;-(6.398+3.5i8iogiox) _ (3 2) 

The signal-to- noise ratio is given, from Eqs. (|2.10|) , and (|2.11|) , by 

/ = 4|^|Vo"'/'/(7)/5o, (3.3) 
where we define the integrals I{q) by 

I{q)= ---dx, (3.4) 

where Xmin = /min/ /o and a;max = /max/ /o, corresponding to the minimum and maximum fre- 
quencies over which the detector will integrate . In some calculations, the maximum value of 
/max corresponds to radiation emitted at the innermost stable circular orbit of the black hole, 
with frequency (6^/^7rm)~^, while in others, we consider the effect of terminating observa- 
tions sooner than this final orbit. The frequency /mm corresponds to the gravitational- wave 
frequency observed a time T earlier, where for LISA-type systems, T < one year. Using the 
quadrupole approximation for radiation damping [Eq. (|1.5| ) with 6 = 0], one can relate the 
frequencies of gravitational radiation at the beginning and end of any time interval T by the 
expression 

/ 256 T s/aX-^*/^ 

«* = «/(^i + — j ■ (3-5) 

We will also wish to study the distance to which such sources can be detected, in order 
to assess the likelihood of detection. Combining Eqs. ( |2.3| ) and ( |3.3| ), we can relate the 
source luminosity distance to source masses and the SNR, 

..^/J^(./o)- 




2-45GpcU^) 1,^) {-][ - I (^1 /(7)V^ (3.6) 
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We first consider a IAMq NS in a quasi- circular inspiral into a BH of mass ranging from 
10^ Mq to 2 X 10^ Mq and an integration time on LISA of one year leading up to the innermost 
stable orbit. We assume a SNR of 10. The parameter errors found are shown in Table |. For 
MBH masses less than about 70,000 Mq, the bound on a; exceeds the current solar-system 
value of 3000, and for low mass black holes {10^ Mq), the bound could be 80 times larger. 
We then determine the bounds on uj for the same range of masses, and for orbits terminating 
at the ISCO, but for shorter integration times, such as 1/2 year or 1/4 year. The results are 
plotted in Fig. 1. Also shown in Fig. 1 are the approximate distances to which such systems 
can be observed. As expected, the bounds weaken almost in direct proportion to integration 
time. Figure 2 illustrates the effect of failing to observe the final orbits leading to the ISCO: 
shown for various BH masses are the bounds for one year of integration terminating a given 
number of days before the ISCO. The rapid fall-off of the bound with time before the ISCO 
indicates the importance of seeing the highest frequency, most relativistic orbits. For sources 
at a given redshift (corresponding to a given Dl), the SNR can be calculated from Eq. (^.61), 
and this used to determine the relevant bound on u. The results for Z = 0.01, 0.03 and 0.05 
are plotted in Fig. 3. To convert from Dl to Z, we use the formula for a matter dominated, 
spatially fiat cosmology Dl = {2/Hq){1 + Z){1 - (1 + Z)-^/"^] with a value of the Hubble 
parameter Hq = 50kms~^ Mpc~^. The signal-to-noise ratios for several illustrative cases are 
shown in the figure. 

In principle, white-dwarf inspirals should give stronger bounds on uj because their sen- 
sitivities are much smaller than that of black holes. However, for small enough black hole 
mass, tidal disruption will require terminating the matched filter before the ISCO, with a 
consequent loss of accuracy in bounding u, as indicated in Fig. 2. We take, as a crude 
measure of the termination of integration for a white dwarf inspiral, the Roche radius 
R ^ d{'mBH / t^wdY^^ -I where d is the radius of the white dwarf, or equivalently, the place 
where the gravitational- wave frequency is given by tt/ = itlbh/^'^ = rnwo/d^. Using the 
approximation, for a non-relativistic WD equation of state, ~ 6 x IQi^ MQ{mw d / Mq)~^/^ , 
and substituting into Eq. ( |3.5| ), we find an approximate cut-off time before the ISCO given 
by 

T,.,_„,^4days(^^^j i^-^j jl-^^Y^j [yJ^) |- 

(3.7) 

When T becomes negative (for mBH > 2.2 x W^Mq, for a 1.4M0 WD), the cut-off is the 
ISCO. Using the Fisher information matrix, we calculate the bounds on u for one year of 
integration of white dwarf inspirals down to the appropriate cut-off time, with a SNR of 
10. For niBH > 2.2 x IO^Mq, the bounds are 2.8 times better than for neutron stars of the 
same mass, as expected from the difference in sensitivities, but for tubh < 2.2 x IO^Mq, the 
improvement gradually decreases because of the loss of orbits near the ISCO. For a IO^Mq 
BH, the bound is only 14 percent better than the corresponding neutron-star bound. 
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B. Follow-up Missions to LISA 



Concepts for a future advanced space-based interferometer have been developed, based 
on the LISA project, which, for want of an official name, we will dub "SuperLISA". The 
motivation for this mission is to detect gravitational waves with a peak sensitivity between 
the ranges of LISA and the ground-based detectors, say around 10~^ Hz. Here, it is hoped, 
there will exist a window free of astrophysically generated sources of background gravita- 
tional waves, such as the white dwarf binaries that contribute to noise in LISA below 10~^ 
Hz. With such a window, it is hoped that cosmological backgrounds of gravitational waves 
may be detectable. 

A hypothetical noise curve for this advanced detector is given by 5'„(/) = Sog{x) with 
X = f/fo, where Q 

So = 2.1 X 10"^^Hz^\ 
/o = 10-iHz, 

^(x) = 0.579x~^ + 1 + 0.064x^ (3.8) 

With this sensitivity, SuperLISA should be able to detect inspirals of intermediate mass 
systems (100 — IO^Mq) to cosmological distances with large SNR. Applying the method 
of the previous subsection to the noise curve of Eq. (p.8|), we find that detection of these 
signals from NS-MBH inspirals may allow SuperLISA to place even more dramatic bounds 
on scalar-tensor theory than LISA. 

We assume integration times of one year terminating at the ISCO and a neutron star 
sensitivity of s^s = 0.2, and consider sources at fixed redshifts of Z = 1, 5, and 10. The 
bounds obtainable from a IAMq NS inspiralling to black holes of a few hundreds of solar 
masses exceed tens of millions, and are dramatically better than the solar-system value for 
essentially all BH masses between 100 and IO^Mq (Fig. 4). For sources at redshift 
Z = 1 and Z = 5, the SNR is well above 10. The drop in the bound on u for Mbh < 150Mq 
is caused by the fact that the initial frequency (one year before the ISCO) for such low- 
mass systems is already around 10~^ Hz, where the instrument has its highest sensitivity; 
most of the data in these cases is being taken against a rapidly increasing background of 
instrumental noise (mostly caused by lack of time resolution related to the instrument's arm 
lengths), hence the bound weakens. 

IV. CONCLUSIONS 

We have found that future observations of inspirals of neutron stars into massive black 
holes by space-based laser interferometric detectors such as LISA may place significant 
bounds on the scalar-tensor coupling parameter uj. For inspirals into black holes as low 
as IO^Mq, the bound could be as large as 240,000 (for a SNR of 10), 80 times larger 
than current solar-system bounds. The bound achievable decreases with increasing black 
hole mass, for a given SNR. A follow-up space interferometer, with a sensitive band at 
intermediate frequencies between those of LISA and the ground based interferometers could 
produce bounds in the tens of millions. 
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The bound on is a strong function of the integration time. For integration times shorter 
than the canonical one year, uo decreases almost in direct proportion, and if the incoming 
wave is not integrated up to the last stable orbit, then the possible bound on uj drops off 
sharply. For neutron-star sensitivities smaller than 0.2, all quoted bounds increase by the 
factor [(0.5 - sns)I^-^^] 

Two important issues have not been addressed. 

The first is our restriction to quasi-circular orbits. This approximation is reasonable for 
stellar mass inspirals of compact objects in galaxies, which are expected to be detectable 
by ground-based interferometers. These are systems where gravitational radiation damping 
has had sufficient time to circularize any pre-existing, two-body eccentric orbit. However, in 
dense galactic nuclei or in dense globular clusters containing massive black holes, compact 
objects are expected to be injected into highly eccentric orbits via interactions with a cloud 
of objects surrounding the hole, and to suffer frequent perturbations by these bodies during 
apholion passage. 

Eccentricity in and of itself should not have a strong effect on the bounds we have 
inferred. A preliminary estimate of the effect of small eccentricity on wave detection shows 
a small drop in SNR for an inspiral with a small eccentricity. For eccentricities of e = 0.1 
or less, the drop in SNR be four percent or less if a quasi-circular inspiral template is used 
to match against an eccentric waveform. For eccentricities of 0.2 to 0.3, the drop in SNR 
will increase to at least 10% and may be as large as 35% for some systems. This drop in 
SNR will be accompanied by a similar drop in detection rates and in the bounds on uj. This 
drop is not too large, however, and can be made even less by using an eccentric waveform 
template (with appropriate BD terms) to match the detected gravitational wave. On the 
other hand, frequent, essentially random perturbations of the orbit of the compact object by 
surrounding stars will make parameter estimation very difficult, even in general relativity. 
How important this effect will be, and whether there exist good data analysis techniques to 
handle it, is a subject of active investigation. 

The second issue is event rate for inspirals into the intermediate mass black holes (100 — 
IO^Mq) that give the best bounds on oj. While enough is known and speculated about the 
rate of occurrence of IO^Mq or lO^M© binary BHs in the centers of galaxies to make such 
systems promising sources for LISA, little to nothing is known about black holes of mass 
WMq or less. Until some estimate can be made of how IO^Mq — IQ'^Mq BH are distributed 
in the Universe, it will also be uncertain whether or not neutron stars can be expected to 
exist around these black holes, and how these neutron stars are expected to behave as they 
inspiral. Nevertheless, there has been recent speculation, based on numerical simulations, 
modeling and some observation, that runaway growth of intermediate mass black holes may 
occur in globular clusters or in young compact star clusters near AGNs pl|j2^ . Another 
mechanism that has been discussed is the gravitational collapse of the first stars, assumed 



to be very massive objects 10 — 10 Mq |23 



It is worth pointing out that testing fundamental theory does not require the same event 
rate that a viable gravitational-wave astronomy does. A rare occurrence of a source with just 
the right characteristics can provide a strong test of scalar-tensor gravity. By comparison, 
while 35 years of pulsar astronomy boasts an average pulsar discovery rate of about 3 per 
month (for a total of over 1000 pulsars), it took only one, PSR 1913+16, to test the general 
relativistic quadrupole for gravitational radiation damping. Still, an estimate of the rate 
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of compact body inspiral into intermediate mass black holes would be desirable, in order 
to evaluate how truly feasible such tests of scalar-tensor gravity using space interferometers 
might be. 
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FIGURES 



Bounds on CO using LISA 



Fixed signal-to-noise =10 




3 4 5 

10 10 10 

Black Hole Mass (M^ 

FIG. 1. Possible bounds on lo vs. BH mass: Three different integration times for LISA are 
shown, one, 1/2, and 1/4 year. A 1.4M0 neutron star with sms = 0.2 and a SNR of 10 are assumed. 
Also shown are approximate distances to the sources, along with the current solar-system bound 
on u of 3000. 
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Bounds on CO using LISA 

Effect of termination before ISCO 

' \ ' \ 




End of Observation (Days before ISCO) 



FIG. 2. Effect of endpoint of integration, on bound on uj, in days before the ISCO. Bounds 
achievable on uo drop off rapidly if signal is not integrated up to last stable orbit. 
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FIG. 3. Possible bounds on uj vs. BH mass as a function of source redshift. Numbers above 
curves indicate approximate SNR for detection. 
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FIG. 4. Bounds on lo vs. BH mass for SuperLISA and sources at redshifts oi Z = 1, 5 and 10. 
A \AMq neutron star with sns = 0.2 is assumed. Numbers above curves indicate approximate 
SNR for detection. Drop in sensitivity for Mbh < ISOMq is because GW frequencies are above 
point of SuperLISA's maximum sensitivity. 
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TABLES 



rriBH 


Ate 




AM/M 




Bound 








(Mo) 




A$e 


(%) 


(%) 


on u 








1000 


2.39 


12.3 


.000192 


.0727 


244549 


.885 


-.996 


-.919 


5000 


3.80 


11.7 


.000544 


.0252 


56326 


.970 


-.998 


-.953 


10000 


5.32 


12.9 


.000768 


.0180 


29906 


.978 


-.997 


-.961 


50000 


26.54 


32.0 


.00243 


.0145 


4780 


.989 


-.998 


-.978 


100000 


103.07 


84.8 


.00610 


.0219 


1458 


.993 


-.999 


-.986 



TABLE I. Estimated parameter errors for 1.4 Mq NS-MBH systems: SNR = 10, integration 
time is one year prior to the ISCO, neutron star sensitivity s^s = 0-2. Error in coalescence time 
tc is in seconds; error in coalescence phase is in radians. Errors in chirp mass ^A and reduced mass 
parameter 77 are in percent. Dimensionless correlation coefficients c'^^ ,c'^^ , and c^^ are also given. 
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